%��糣��
% d1_1=0.823 d1_2=0.597 d2_1=0.805 d2_2=0.692 d3=1.689 d4=1.678
clc;
clear;
result1_1 = readmatrix('sam1-1-results.csv');    
result1_2 = readmatrix('sam1-2-results.csv');
result2_1 = readmatrix('sam2-1-results.csv');
result2_2 = readmatrix('sam2-2-results.csv');
result3 = readmatrix('sam3-results.csv');
result4 = readmatrix('sam4-results.csv');
f = result1_1(:,1);
e1_1 = result1_1(:,4);
e1_2 = result1_2(:,4);
e2_1 = result2_1(:,4);
e2_2 = result2_2(:,4);
e3 = result3(:,4);
e4 = result4(:,4);

%%
%�ֿ���
figure('Renderer', 'painters', 'Position', [500 0 420*2 560*3],'color',[1 1 1])
subplot(321);plot(f,e1_1,'d-','color',[107, 45, 35]/255,'LineWidth',1.5,'MarkerSize',6,'MarkerIndices',1:4:length(e1_1),'MarkerEdgeColor',[107, 45, 35]/255);
axis([0.2 2 2.6 3.1])
set(gca,'XTick',0.2:0.3:2);
set(gca,'FontSize', 14,'FontWeight','bold');
legend('d=0.823mm','Location','northeast')
box on;

subplot(322);plot(f,e1_2,'h-','color',[242, 174, 37]/255,'LineWidth',1.5,'MarkerSize',6,'MarkerIndices',1:4:length(e1_2),'MarkerEdgeColor',[242, 174, 37]/255);
axis([0.2 2 2.6 3.1])
set(gca,'XTick',0.2:0.3:2);
set(gca,'FontSize', 14,'FontWeight','bold');
legend('d=0.597mm','Location','northeast')
box on;

subplot(323);plot(f,e2_1,'o-','color',[202, 102, 26]/255,'LineWidth',1.5,'MarkerSize',6,'MarkerIndices',1:4:length(e2_1),'MarkerEdgeColor',[202, 102, 26]/255);
axis([0.2 2 2.6 3.1])
set(gca,'XTick',0.2:0.3:2);
set(gca,'FontSize', 14,'FontWeight','bold');
ylabel('Dielectric constant','FontSize', 20,'FontWeight','bold')
legend('d=0.805mm','Location','northeast')
box on;

subplot(324);plot(f,e2_2,'d-','color',[173, 58, 28]/255,'LineWidth',1.5,'MarkerSize',6,'MarkerIndices',1:4:length(e2_2),'MarkerEdgeColor',[173, 58, 28]/255);
axis([0.2 2 2.6 3.1])
set(gca,'XTick',0.2:0.3:2);
set(gca,'FontSize', 14,'FontWeight','bold');
legend('d=0.692mm','Location','northeast')
box on;

subplot(325);plot(f,e3,'h-','color',[147, 16, 14]/255,'LineWidth',1.5,'MarkerSize',6,'MarkerIndices',1:4:length(e3),'MarkerEdgeColor',[147, 16, 14]/255);
axis([0.2 2 2.6 3.1])
set(gca,'XTick',0.2:0.3:2);
set(gca,'FontSize', 14,'FontWeight','bold');
xlabel('Frequency(THz)','FontSize', 18,'FontWeight','bold')
legend('d=1.689mm','Location','northeast')
box on;

subplot(326);plot(f,e4,'o-','color',[115, 7, 7]/255,'LineWidth',1.5,'MarkerSize',6,'MarkerIndices',1:4:length(e4),'MarkerEdgeColor',[115, 7, 7]/255);
axis([0.2 2 2.6 3.1])
set(gca,'XTick',0.2:0.3:2);
set(gca,'FontSize', 14,'FontWeight','bold');
xlabel('Frequency(THz)','FontSize', 18,'FontWeight','bold')
legend('d=1.678mm','Location','northeast')
box on;